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It is possible using current high intensity laser facilities to reach the quantum radiation reaction 
regime for energetic electrons. An experiment using a wakefield accelerator to drive GeV electrons 
into a counterpropagating laser pulse would demonstrate the increase in the yield of high energy 
photons caused by the stochastic nature of quantum synchrotron emission: we show that a beam 
of 10 9 1 GeV electrons colliding with a 30 fs laser pulse of intensity 10 22 Wcm -2 will emit 6300 
photons with energy greater than 700 MeV, 60 x the number predicted by classical theory. 



FIG. 1. (Color online). Diagram of an experimental geometry 
that could demonstrate quantum radiation reaction and pair 
production. The GeV electrons decelerate in large fields at 
the laser focus, producing gamma rays that pass through a 
hole in the f/2 optic. 


Now that the power focussed in laser facilities exceeds 
1 PW, electron dynamics enters a regime between classi¬ 
cal and QED physics. In this letter we will consider an 
experimental setup where this transition can be explored: 
the collision of a GeV electron beam and a high-intensity 
laser pulse operating close to the current intensity fron¬ 
tier ( 10 22 Wcm -2 ), as shown in Fig. [lj Electron dynam¬ 
ics in this regime is dominated by radiation reaction [l|, Q, 
which will manifest itself in the efficient conversion of 
electron energy to hard gamma rays. We show that with 
current high-intensity laser facilities, it is possible to pro¬ 
vide the first demonstration of the probabilistic nature of 
radiation reaction in the multi-photon, strong-field QED 
regime. The stochastic nature of this process increases 
its efficiency by over an order of magnitude compared to 
classical radiation reaction, and provides the most sensi¬ 
tive diagnostic of quantum synchrotron emission. 

This experimental setup will provide a means of testing 
the fundamental physics underlying more exotic phenom¬ 
ena, such as pair cascades 0 0, which occur in pulsar 
magnetospheres [h. 6] and are predicted to become signif¬ 
icant above a threshold between 10 23 and 10 24 Wcm -2 . 
In these cases the plasma dynamics is strongly affected 
by the generation of macroscopic electron-positron pair 


plasmas, for example at the laser focus in laser-solid tar¬ 
get interactions [zj. 

We parameterise the importance of strong-field QED 
effects with the quantity ij = \F llu p u \/mcEs c h — 7|E_l + 
v x B|/Esch 0, for an electron with four-momentum 
pH _ ^ m ( C5 v ) i n an electromagnetic field F^ v with mag¬ 
netic component B and component of the electric field 
perpendicular to the electron direction of propagation 
E ± . r] is the ratio of the electric field in the electron 
rest frame to Es c h = m 2 c 3 /eh, the characteristic field of 
QED [ 9 ]. In the setup we consider, the probabilistic na¬ 
ture of photon emission is evident for 77 ~ 0 . 1 . The QED- 
dominated regime is reached when r] > 1 ; at 77 = 90 , a 
50 GeV electron would be capable of producing multiple 
pairs on passage through a petawatt laser pulse 0- 

In a classical description, the electron travels on its 
worldline radiating continuously as it accelerates in the 
laser fields. The typical energy of a single emitted pho¬ 
ton is 0.44?77mc 2 [S|; therefore as r] approaches 1 , a single 
photon becomes capable of carrying off a significant frac¬ 
tion of the electron’s energy and the recoil of that emis¬ 
sion must be taken into account. However, the quantum 
description of radiation reaction differs from classical the¬ 
ory in two ways. 

Firstly, quantum corrections to the radiated spectrum 
cut off the tail of photons with energies greater than that 
of the electron [ljJ and include spin-flip transitions, the 
latter markedly increasing the probability of radiating 
photons with energies comparable to that of the electron 
[ 12 J. These modifications means that the total radiated 
power is smaller than the equivalent classical power by a 
factor g(rj) G ( 0 , 1 ) jH. fl^. 

The second and more significant effect for us is that 
the process of photon emission is stochastic 0, [HI; thus 
the electron has only a probability to emit a gamma-ray 
photon of given energy. This gives rise to a phenomenon 
called ‘straggling’ [l 5 |, ll6(| , where the electron may propa¬ 
gate a significant distance through the strong laser fields 
without radiating. Since the laser pulse will have a spa¬ 
tial intensity profile, it is possible for the electron to reach 
the region of highest intensity at the centre having lost 
much less energy than a classical electron. The 77 of an 
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electron that has straggled in this way will be boosted 
above that which could be reached classically. As the 
tail of the photon spectrum increases non-linearly with 77 , 
straggling enhances the yield of hard gamma-rays. Fur¬ 
thermore, these energetic photons can decay in the strong 
electromagnetic fields to produce electron-positron pairs 
by the Breit-Wheeler process [s]. Pairs can also be pro¬ 
duced directly by the electron in the trident process. 

These strong-field QED effects have been investigated 
experimentally with 100 GeV electrons incident on crys¬ 
tals at the CERN Super Proton Synchrotron 0 and 
in the collision of a 50 GeV electron beam and a 
10 18 Wcm -2 laser pulse at the SLAC facility jTsj. In the 
former, the reduction in radiated power g(rj) was mea¬ 
sured and found to agree well with theoretical predic¬ 
tions; in the latter, pair production by gamma rays cre¬ 
ated by inverse Compton scattering was observed. How¬ 
ever, by using a more intense laser, we can reach the 
same 77 with electrons of much lower energy. As it has 
been shown that PW lasers in long focus can generate 
GeV electron beams by wakefield acceleration 19l-[21| . it 
is realistic to consider the experiment shown in Fig. HJ 
where another laser pulse, in tight focus, provides the 
high-intensity target for the wakefield-accelerated elec¬ 
tron beam. 

We have developed a Monte-Carlo algorithm to simu¬ 
late the interaction of an intense laser pulse with an en¬ 
ergetic electron beam, following §2 of [15] and §3 of [22]. 
The three processes included are: photon and (trident) 
pair production by an electron, and (Breit-Wheeler) pair 
production by a photon. 

The rates of these three processes are calculated in the 
Furry picture of QED [23], where the electron interacts 
with both an external, unquantised electromagnetic field 
and a fluctuating component of the same [24. 25]. In be¬ 
tween interactions with the latter (i.e. photon emission 
and absorption) it propagates classically [26]. The pho¬ 
ton formation length is smaller than the laser wavelength 
by a factor of the laser strength parameter ao (27] ; since 
a 0 = [/L(A//im) 2 /1.37 x 10 18 Wcm -2 ] 1 / 2 >> 1, where Jl 
and A are the laser intensity and wavelength respectively, 
we treat the emission process as pointlike and instanta¬ 
neous. As the fields are quasi-static over the emission 
process, these rates may be calculated in an equivalent 
system of fields with the same instantaneous value of 
77 , such as a static magnetic field (in the limit B 0, 
7 —>> 00 ) or a plane EM wave (in the limit of zero fre¬ 
quency) [ItJ. 

Using the static magnetic field approach, the spectrum 
of emitted photons is determined by the quantum syn¬ 
chrotron function F(t 7 ,x), where x = ( h/2mc 2 )\uj'E± + 
c 2 k x B|/Esch is the counterpart of the electron parame¬ 
ter 77 for a photon of frequency uj and wavevector k. The 
rate of emission is given in Erber 0 and Baier et al. 
0 . In the plane wave approach, this same process is 
called non-linear inverse Compton scattering 26, 28, 29j . 


We further assume that the electrons may be treated 
independently as they are highly energetic, with 7 
ao ^ 1. The rigid-beam approximation is valid as the 
transverse momentum gained from the laser is small, and 
in the frame co-moving with a nearly mono-energetic 
beam, the energy of interaction between particles will 
also be small. Therefore the electron motion, initially 
antiparallel to the optical axis, is one-dimensional and 
constant between discrete photon emission events. 

This algorithm is implemented thus: the code initially 
assigns to each electron a pseudorandom ‘final’ optical 
depth against photon emission r p h and trident pair pro¬ 
duction 7 ”t r i. As they propagate through the laser pulse, it 
integrates their differential optical depths against trident 
pair production dr tr i/dt and photon emission dr p h/dt = 
(y/3a/27rrc)rjh(rj)/ r y, where a is the fine structure con¬ 
stant, re the Compton time and h(rj) = J ^ 2 dxF(r], x)/x 

HEE9. 

Emission occurs when the electron reaches that ‘final’ 
optical depth [uj|. Following that, a new final depth is as¬ 
signed. The photon energy is found by pseudorandomly 
sampling the quantum synchrotron distribution d, 10 ] 
(J 0 x dxF( 77 , x)/x)/h(rj). The electron energy, constant 
between emissions, is then reduced by the energy of the 
emitted photon. In the static magnetic field approach, 
this follows from energy conservation. Since synchrotron 
radiation does not conserve the component of momentum 
normal to the magnetic field [30|, when implemented in 
the code, this leads to a fractional error A 7/7 oc 1 / 7 em it 
which is negligible for all emission events [4|, |2| . 

Energy loss to trident pair production is neglected, as 
the rate is calculated in the Weizsacker-Williams approx¬ 
imation, which treats the exchanged virtual photon as 
real (i.e. the momentum transfer from electron to pair 
q 2 0). Breit-Wheeler pair production is modelled by 
assigning to each photon at its creation a final optical 
depth tbw- The rate of pair production [b3] is integrated 
along the photon trajectory to determine whether it de¬ 
cays to a pair. 

We will compare our fully stochastic model of radia¬ 
tion reaction with a semi-classical model, which we call 
‘continuous radiation reaction’. In this model the elec¬ 
tron loses energy continuously according to a damping 
term given by the Landau-Lifshitz force |30| modified to 
include the quantum correction to the radiated power: 
dj/dt = (2a/3)(r] 2 /rc)g(r]). For consistency, photon 
spectra are obtained by sampling at each timestep the 
quantum synchrotron distribution. Thomas et al. [3l[ 
use a similar model to simulate laser-electron collisions 
for electrons with 70 = 400, where 77 is small enough that 
probabilistic effects can be neglected. 

The parameters of the simulation are as follows: the 
laser pulse has wavelength A = 1 fi m, is linearly polarised 
and has Gaussian temporal profile with a full width at 
half-maximum (FWHM) of 30 fs. The electrons have ini- 
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FIG. 2 . (Color online). The number of photons with en- -S 

ergy > ftw for a 1 GeV electron incident on a linearly po- ^ 

larised, plane-wave laser pulse with Gaussian temporal profile 
(FWHM 30 fs) and peak intensity (a) 10 21 (b) 10 22 Wcm -2 . 

Inset (top) The fractional increase, due to straggling, in the 
number of photons with given minimum energy. Inset (bot¬ 
tom) The maximum 77 experienced by a (c) continuously and 
(d) discontinuously radiating GeV electron incident on the 
same pulse with given intensity. The increase in going from 
(c) to (d) is responsible for the hardening of the photon spec¬ 
trum. 

tial gamma factor 70 = 2 x 10 3 and propagate along the 
optical axis antiparallel to the laser pulse; this ensures 
that the and v x B terms in the definition of 77 are 
additive, i.e. 7 ~ 2^/Ei,/Es c h- Each simulation follows 
at least 10 7 macroelectrons, their trajectories discretised 
into intervals of At = 10 _ 17 s. The laser fields are calcu¬ 
lated at each timestep using the classical solution for a 
Gaussian focussed beam. 

If the laser pulse is a plane wave, the number of pho¬ 
tons per electron with energy greater than ftcj, and the 
increase in the same due to straggling, are given in Fig. 

[2j The boost to the electron’s maximum 77 as a result of 
straggling is also shown: this boost is responsible for the 
hardening of the photon spectrum. 

In reality both the laser and electron beam will have 
temporal and spatial structure. As long as enough of the 
beam is incident on the laser focus, there will still be 
strong evidence of probabilistic photon emission. Con¬ 
sider a beam of 1 GeV electrons propagating antiparallel 
to the laser but uniformly distributed around the optical 
axis in a disk of radius 10 fi m. The laser pulse, still lin¬ 
early polarised, is focussed with Gaussian temporal and 
radial profiles (FWHM 30 fs and waist size 2.2 /1 m). If it 
has peak intensity of 10 22 Wcm -2 , the number of photons 
per electron with energy greater than 500 (700) MeV this 
interaction produces will be 7.40 x 10 -4 (1.32 x 10 -5 ), 

4.3 x (160 x) greater than that which would be radiated 
classically. 

It is then possible to distinguish discontinuous radi¬ 
ation reaction by measuring the high-energy tail of the 
gamma-ray spectrum, as the yield of these photons is 



1.0 1.5 2.0 3.0 5.0 7.0 10.0 15.0 20.0 30.0 

Peak intensity (10 21 Wcm -2 ), cont. rad. 

FIG. 3. (Color online). Indicated along the coloured lines are 
the numbers of photons per electron with energy greater than 
100 MeV (red), 200 MeV (green), 400 MeV (blue) and 500 
MeV (purple) if the electrons radiate discontinuously (vertical 
axis) or continuously (horizontal axis). The interaction here is 
between a focussed laser pulse of given peak intensity (FWHM 
30 fs, waist 2.2 fi m) and a beam of 1 GeV electrons (radius 
10 fim around the optical axis). 


most enhanced by straggling. Fig. [3] shows that if the 
electrons radiate semi-classically, rather than stochasti¬ 
cally, the laser pulse must be much more intense for the 
interaction to produce the same number of high-energy 
photons. This increase in intensity is necessary because 
a stochastically radiating electrons can reach a higher 77 . 

Assessing the nature of the emitted radiation in a laser- 
electron beam experiment can then done by fitting the 
observed spectrum to the points in Fig. 03 For ex¬ 
ample, the detection of 9.0 x 10 -2 photons per elec¬ 
tron with energy greater than 100 MeV would be con¬ 
sistent with either a continuously or discontinuously ra¬ 
diating beam colliding with a laser pulse of peak inten¬ 
sity 5 x 10 21 Wcm -2 . In contrast, the simultaneous de¬ 
tection of 1.9 x 10 -3 photons per electron with energy 
greater than 400 MeV would only be consistent with dis¬ 
continuous radiation reaction, as a continuously radiat¬ 
ing beam would have to collide with a pulse of peak in¬ 
tensity 1.8 x 10 22 Wcm -2 that number of gamma rays. 
If the number of photons with energy greater than 500 
MeV is measured instead, that equivalent intensity rises 
to 5.2 x 10 22 Wcm -2 . 

A typical laser wakefield will accelerate a bunch of 10 9 
electrons, which will produce a beam of 4.1 x 10 5 photons 
with fajj > 500 MeV, collimated within an angle of 30 
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FIG. 4. (Color online). The number of Breit-Wheeler (red) 
and trident (blue) electron-positron pairs produced per 1 GeV 
electron colliding with a linearly polarised plane-wave laser 
pulse with Gaussian temporal profile (FWHM 30 fs) and given 
intensity. Inset The fractional increase in the number of pairs 
for the same processes produced by treating radiation reaction 
probabilistically. 


mrad. When shielded to block lower energy gamma rays, 
this will be sufficient energy deposition in a calorimeter to 
provide a good signal of discontinuous radiation reaction. 

Our code also models pair production by both the 
gamma rays and the electrons. Shown in Fig. 0 ] are 
the number of pairs produced per electron by a monoen- 
ergetic 1 GeV beam incident on a plane-wave laser pulse 
with given peak intensity. 

Pair production will form a key part of the new physics 
probed by future high-intensity laser experiments. How¬ 
ever, pairs do not become energetically significant in the 
experimental setup we describe here because they are not 
accelerated by the laser fields; thus their probability of 
emitting a photon or creating an additional pair before 
exiting the laser pulse is very small. Converting laser en¬ 
ergy to pairs requires much more intense, counter propa¬ 
gating laser beams that can accelerate the electron up to 
high energy between emission events (i. 15. E3- 

The large increase in the number of Breit-Wheeler 
pairs shown in the main part of Fig. 0 ] arises because the 
rate drew/d t is highly non-linear in x an d has maximum 
rate of growth for x = 0.62. Reaching this switch-on with 
current high-intensity lasers is strongly dependent on the 
yield of close to GeV photons. As straggling enhances 
this yield, it also increases the yield of pairs (see inset 
of Fig. 0 ]). The number of trident pairs is similarly in¬ 
creased, as straggling allows more high-energy electrons 
to penetrate to the region of highest field intensity. 

Although potentially more difficult to measure, the 
electron energy distribution can provide evidence of 
quantum radiation reaction [3^]. In Fig. [5j we show how 
the initally monoenergetic beam acquires a spread be¬ 
cause the electrons counterpropagate with a range of dis- 
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FIG. 5. (Color online). The final state energy distribution 
dN/dE e of a 10 /un radius beam of 1 GeV electrons incident 
on a pulse with given peak intensity. White, dashed lines 
are contours of constant dN/dE e . The red line is the lowest 
energy that can be reached with classical radiation. 


placements from the optical axis, experiencing different 
laser intensities. The energy loss of electrons subject to 
continuous radiation reaction is determined only by laser 
intensity. Thus there is a non-zero lower bound to the 
classical final state energy distribution, corresponding to 
those electrons which have passed through the region of 
highest laser intensity. 

Discontinuously radiating electrons lose energy prob¬ 
abilistically, so some electrons will straggle and lose 
more energy than possible classically. At 10 22 Wcm -2 , 
6.1 x 10 -3 of the beam electrons will experience an en¬ 
ergy loss exceeding 863 MeV, the maximum possible clas¬ 
sically. It may be possible to distinguish these electrons 
if the initial energy of the beam is very well-characterised 
and devoid of a low energy tail. This could be accom¬ 
plished with a particle accelerator 0 , or with magnetic 
filtering of the electron beam produced by a laser-driven 
wakefield. 

In this letter we have considered the effects of including 
a fully stochastic model of radiation reaction on the mo¬ 
tion of an energetic electron beam incident on an intense 
laser pulse. Under the chosen conditions, the parameter 
Tj = 27FV/£’sch ~ 0.1 and so QED effects must be in¬ 
cluded. We find that electron motion is dominated by 
quantum radiation reaction, enhancing the yield of pho¬ 
tons with Huj > 700 MeV by a factor of 160; thus the col¬ 
lision of a GeV electron beam with a petawatt laser can 
provide evidence of both stochastic gamma ray and pair 
production. The observation of 700 MeV photons will be 
an unambiguous signature that GeV electrons have been 
incident on, and straggled through, the region of high¬ 
est intensity at the laser focus. The measurement of in¬ 
creased energy spreading of the electron beam is likely to 
be more ambiguous in a realistic experiment where laser 
parameters vary from shot to shot. Experimental valida¬ 
tion of the probabilistic nature of QED-plasma models 
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will underpin the simulation and design of the next gen¬ 
eration of higher energy laser-plasma interactions. 
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